Griffiths phase in the thermal quantum Hall effect 
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Two dimensional disordered superconductors with broken spin-rotation and time-reversal invari- 
ance, e.g. with p x + ip y pairing, can exhibit plateaus in the thermal Hall coefficient (the thermal 
quantum Hall effect) . Our numerical simulations show that the Hall insulating regions of the phase 
diagram can support a sub-phase where the quasiparticle density of states is divergent at zero en- 
ergy, p(E) ~ li?! 1 / 2-1 , with a non-universal exponent z > 1, due to the effects of rare configurations 
of disorder ("Griffiths phase"). 

PACS numbers: 73.43.-f; 73.20.Fz; 75.10.Nr 
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The integer quantum Hall effect (IQHE) is observed 
in two-dimensional (2d) electron gases in high magnetic 
fields. Its striking manifestation is the quantization of 

2 

the Hall conductance o~ X y in units of t- . Analogs of the 
IQHE are also known in systems with broken time re- 
versal symmetry in which quasiparticle charge conserva- 
tion is violated because of the presence of a charge con- 
densate, while quasiparticle spin and energy remain con- 
served. In such a situation the Hall coefficients for spin 
or heat transport can exhibit plateaus at integer values 
in units of h/2-ir Q or Trk 2 B /3H 0,0, respectively. These 
novel types of quantum Hall effects may be encountered 
in superconducting systems with unconventional (non- 
s-wave) pairing. There has been a great deal of recent 
interest in these spin and thermal quantum Hall effects 
(SQHE and TQHE) [JSSSEllHH- In P artic " 
ular, the TQHE 9] in systems without spin conserva- 
tion, which belongs to class D in the Altland-Zirnbauer 
scheme , displays many peculiar features that are not 
found in the IQHE or SQHE, including a metallic phase 
with a logarithmically diverging density of states (DoS) 
0, [n| . More generally, behavior of the DoS in uncon- 
ventional disordered superconductors has attracted a lot 
of research interest recently 12] . A wealth of different 
types of low-energy behavior of the DoS was identified, 
depending on the symmetry of the system and on the 
type of disorder. 

In a recent study of the Random Bond Ising model 
(RBIM) (which falls in the same universality class as the 
TQHE systems), Gruzberg et al. conjectured that 
such systems may have Hall insulating phases support- 
ing a region where the DoS exhibits a power law singu- 
larity p(E) oc |£'|~ 1+1 / z with a nonuniversal exponent z. 
Strong-disorder renormalization group (RG) calculations 
of Motrunich et al. [l4| have also demonstrated the pres- 
ence of such low E divergence in the DoS (with z > 1) in 
quasi- Id superconductors with broken time reversal and 
spin rotation invariance. This has also been confirmed by 
calculations using a Fokker-Planck approach [J^. The 
strong-disorder RG calculations make it apparent that 



the DoS divergence originates in this and in other Id 
models |l6| from Griffiths effects involving exponen- 
tially rare length I regions (occuring with probability 
~ e~ c ) of the sample in 'nearby' phases. An exponen- 



tially weak coupling, ~ e , of defects at the boundaries 
of such rare regions produces a power-law tail in the low 
energy DoS: 



P(E) 



Jail 5(E -e- cl ) e~ c ' 1 oc 



(1) 



[z=c/d and c, c' are non universal parameters). 
Motrunich et al. further speculated [Til EH that a simi- 
lar mechanism driving a divergent DoS may be operative 
also in 2d superconductors. 

The presence of such a regime with divergent DoS is 
in sharp contrast to the situation in standard Wigner- 
Dyson ensembles, where the DoS is never divergent. In 
2d power law divergences generically modified by loga- 
rithmic corrections have been known to exist in the chi- 
ral symmetry classes [l8|, [19J, |2fj, |21| , where they appear 
in critical phases without exponential localization of the 
wavefunctions. Here, we report an observation of an in- 
sulating Griffiths phase with power-law singularities in 
the DoS of a model (RBIM) in the TQHE universality 
class. We further demonstrate that the divergent DoS is 
governed by a mechanism involving the presence of rare 
long string defects jl^ . Our numerical results thus show 
that TQHE - in contrast to its close relatives, the IQHE 
and the SQHE - does exhibit the conjectured Griffiths 
phase (HQ. 

For our numerical simulations we choose a particular 
network realization of the ± J RBIM. As we discuss in the 
end of the paper, our results are however not restricted 
to this model but rather are generally valid for TQHE 
systems. The RBIM network realization is defined on a 
square lattice and has been described in detail before [jj . 
Particles traversing the links of the network change di- 
rection and turn left or right with probability amplitudes 
±cosa or ±sina whenever arriving at a node. (The re- 
spective signs are fixed by the condition of unitarity of the 



2 



LT 



N / 



tie »^>I«»I < 

Griffith 

phase*' ferro 




para 



HT 



FIG. 1: Phase diagram of TQHE in RBIM-network model 
representation, disorder p vs. inter plaquette coupling 
strength (temperature) ti. Fixed points: high temperature 
(HT), low temperature (LT), clean Ising (I), multicritical 
point (N), percolation (P c ). Phase boundary (solid) sepa- 
rates two Hall insulators (called ferro and para using RBIM 
terminology) . Shaded region indicates Grifhths phase with di- 
verging single particle density of states, p(E). Line segments 
represent scans for p(E) displayed in Figs. [5] (x) and[3](*). 
In Fig. El scan is along phase boundary from N to P c . 



corresponding time evolution operator.) Disorder in the 
form of "vortex pairs" is introduced by flipping the sign of 
both cos a amplitudes of the node that is dressed by the 
vortex pair. The system is characterized by two param- 
eters: the inter-plaquette coupling strength, ti = sin 2 a 
related to the temperature T and coupling constant J 
of the Ising model as ti = cosh -2 2 J/T, and the con- 
centration of defects, p. Note that although the inter- 
plaquette hopping amplitudes are very different from the 
intra-plaquette amplitudes, the underlying spin model 
has no explicit dimerization in its couplings. Fig^shows 
the phase diagram of the model. It contains two local- 
ized phases corresponding to the ferro- and paramagnetic 
phases of the RBIM. The precise location of the phase 
transition line (the thermal quantum Hall transition in 
the fermionic formulation) has been studied in Ref. |25| . 
We also show in Fig. ^the Nishimori line (NL), where 
the RBIM possesses some additional symmetry. 

We construct the time evolution operator U for a single 
time step for networks with periodic boundary conditions 
and use standard sparse matrix packages to extract a 
number of eigenvalues closest to zero |26j . Typically 8 
to 16 eigenvalues per system are calculated. Since U 
is unitary and the network is class D, the eigenvalues 
come in pairs and are of the form exp(±iE). We obtain 
the density of (pseudo)-energies, p(E), from an average 
over an ensemble of disorder realizations which contains 
typically of the order of 10 4 members. 

We start the presentation of our numerical results with 
a scan along the line p = 0.08, as shown in Fig. ^ The 
scan begins in the bigh-T (paramagnetic) phase where 
the impact of a small concentration of vortex disorder is 
relatively weak. Here, one expects the energy gap that 
exists in the clean limit (p=0) to dominate because the 
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FIG. 2: DoS near (top) and on (bottom) of phase boundary at 
p=0.08 (see Fig. 0. Top: Bare gap with exponential smearing 
by Lifshitz mechanism. Numerical fit of gap region is compat- 
ible with In pocE~ K , k w 2. Bottom: Pseudo-gap in DoS fitted 
to the analytical prediction p(E)Ri(g M /2n 2 )\E\ In \Ei/E\; fit: 
3M-21.51 and £1 =0.262. 

particles have an appreciable tendency to wind about 
those plaquettes of the sublattice that are defect free. 
A non-zero DoS near zero energy is only due to very rare 
configurations of vortices giving rise to a "Lifshitz tail " 
[13, In p(E) ~ E^ 2 , see top panel of Fig.^that turn the 
true gap of the clean system into a "pseudo-gap" . 

At lower values of the interplaquette coupling, t\ w 
0.29, i.e. right on the phase boundary, the shape of the 
DoS is more complicated. In the absence of disorder, the 
DoS is linear, p(E) oc \E\, at the Ising critical point. In 
a disordered system, the theory of Dirac fermions with 
a random mass predicts a logarithmic correction to the 
low-energy behavior, p(E) cx \E\ In [ll|. Indeed our 

data for not too high E are consistent with this form, see 
bottom panel of Fig. (at higher energies model-specific 
band structure effects become important). Small devia- 
tions from this formula at very low energies can be at- 
tributed to finite-size effects, which are particularly pro- 
nounced at criticality. 

If we diminish the inter-plaquette coupling strength so 
as to get close to the NL (still keeping p = 0.08 constant), 
we witness a qualitative change in the low-energy DoS: 
the pseudo-gap vanishes and the DoS starts to increase 
with lowering E, apparently diverging as E — > 0. This 
is demonstrated in Fig. where the evolution of DoS 
is shown as ti is lowered from 0.13 down to 0.014. For 
large enough inter-plaquette coupling one still notices the 
pseudo-gap behavior. The latter is already less evident 
in the data that corresponds to the coupling strength 
ti = 0.07. It finally gives way to a strongly increasing 
behavior once the coupling falls below ti w 0.05. 

Having demonstrated the divergence of the DoS in the 
low-temperature part of the ferromagnetic phase, we turn 
to a discussion of the corresponding mechanism. We will 
show that the DoS singularity is governed by rare disor- 
der configurations containing long chains (strings) of de- 
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FIG. 3: DoS at disorder value p = 0.08 deep in the ferromag- 
netic phase. As coupling between plaquettes is progressively 
reduced pseudo-gap gives way to power law divergency. Plots 
show averages for two system sizes, 128 (256) with 20000 
(10000) realizations. Data comprises 8 lowest lying eigen- 
states; t\ = 0.029 corresponds to location on Nishimori line. 

fccts. Since each defect contains a pair of vortices, such 
a string of a length £ > 1 (formed with a small but non- 
zero probability ~ p e ) produces two vortices with a large 
separation £. At zero plaquette coupling, ti = 0, every 
vortex contributes a zero mode (or "mid-gap state") to 
the spectrum which therefore has a macroscopic degen- 
eracy pL 2 . Upon switching on the coupling ti, the zero 
modes couple and the degeneracy is lifted. For a small 
ti , the energy splitting E thus produced will be exponen- 
tially small in I j^. For a clean system (which contains 
only two vortices at the ends of the string but no other 
defects) this follows immediately from the perturbation 
theory, yielding E oc t l J 2 = e~ ct . Clearly, the same result 
holds for those rare configurations of the disordered sys- 
tem for which the region of a size ~ £ around the vortex 
pair is free of other vortices. However, the probability 
for this to happen scales with I as (1 — p) , where d = 2 
is the spatial dimensionality. In the Id case this factor 
is harmless, yielding the probability for formation of an 
isolated string oc e~ c 1 with d = — lnp(l — p) 2 and Eq. Q] 
applies. On the other hand, in d = 2 the above factor 
seems to spoil the reasoning leading to Eq. QJ. So, if 
the requirement of "isolation" of the vortex pair would 
be a necessary condition for the appearance of exponen- 
tially small eigenvalues, e ~ e~ ci , the divergent behavior 
of the DoS would be converted to a pseudo-gap at lowest 
energies due to the (1 — p) factor. 

We now argue that this requirement is not necessary, 
and Eq. (|TJ is applicable in 2d as well. Indeed, at small 
p almost all randomly located defects will be isolated. 
Correspondingly, ~ p is the concentration of vortex pairs 

1/2 

with separation unity and an energy splitting ~ t, ' . 
There will be a much smaller concentration (p 2 ) of de- 
fects forming chains of length two, yielding the energy 



FIG. 4: Distribution of logarithm of lowest eigenvalue for 
string of defects of length I = 8 -r 56 inserted in disordered 
system in ferromagnetic phase (p — 0.08, ti — 0.029, L = 
128). Inset: typical lowest eigenvalue as a function of I. 

splitting ~ ti, and so forth, the energy splitting as- 
sociated with these short (compared to t) vortices will 
be relatively big, and will only weakly affect the energy 
splitting t 1 -! 2 of the ^-vortex. This arg ument can be con- 
verted to a RG procedure, as in [14i l21| . The vortices 
can be sorted into vortex pairs by combining those with 
the shortest distance first, the second shortest distance 
next, etc. Since the energy splitting associated with the 
coupling of two vortices is exponential in their distance, 
vortex pairs differing in this distance will strongly differ 
in their splitting as well. The hierarchical order allows 
for a RG procedure that successively eliminates vortex 
pairs with an energy splitting that is large as compared 
to the pairs remaining on the lattice. Clearly, for suffi- 
ciently small ti and p (the boundary of the Griffiths phase 
will be discussed below) the eigenvalues associated with 
long strings, i.e. vortex pairs with a large separation, are 
(largely) insensitive to the presence of other vortex pairs. 
In other words, a string of length I does not have to be 
isolated in order to contribute an eigenvalue ~ e~ cl . The 
only role of other (randomly distributed) vortices will be 
in somewhat modifying the parameter c as compared to 
its value c = — g lnt; for an isolated string. 

Since the validity of the above RG-type argument is 
crucial for the proposed string mechanism, we have per- 
formed a numerical test for its verification. We have in- 
serted a string of the length i (ranging from 8 to 56) into 
a disordered network with t; = 0.029 and p — 0.08. Fig. 
0] displays the evolution of the distribution of the lowest 
eigenvalue with increasing length of the string. The re- 
sults fully confirm our scenario: a) the extra string gives 
rise to an additional, isolated eigenvalue; b) this eigen- 
value moves with increasing / according to the expected 
exponential law, (— lne) oc £, see inset of Fig. 01 (the 
rest of the spectrum does not move). We have performed 
the same calculation in the paramagnetic phase (p = 0.2, 
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that occurs in the TQHE characterized by a divergent 
density of states. Our numerical work is based on the net- 
work representation of the ± J random bond Ising model. 
Preliminary results for the more generic Cho-Fisher |2flj 
model confirm that a similar phase exists also there, so 
that the Griffiths phase is a ubiquitous companion of the 
TQHE. 
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FIG. 5: Scan for DoS along phase boundary starting from 
Nishimori point (N) terminating at P c . DoS shows power- 
law divergence, p(E) ~ |_E| 1/ ' Z ~ 1 , with non-universal ex- 
ponent z. At low energies random-matrix oscillations in- 
dicate existence of delocalized (critical) states. Plots with 
t t = 0.058,0.044,0.029 (from top to bottom); each for three 
values of system size: L — 64, 128, 256. 

ti = 0.029) and found that the distribution does not move 
with I: a long string of defects in a random environment 
does not generate an exponentially small eigenvalue. 

Therefore, the analytical arguments together with the 
numerical results, Figs. and 0] demonstrate the exis- 
tence of a Griffiths phase where string formation gener- 
ates a divergent DoS, Eq. with z > 1. This phase 
is shaded in the phase diagram, Fig. ^ Its boundary 
on the side of larger ti is given by the condition z = 1 
and, according to our numerical results, is at least close 
to the NL. The asymptotic behavior of the boundary at 
low temperature and weak disorder can be easily found 
from our analytical considerations. Indeed, in this limit 
c ~ — | In i<> ~ 2J/T, and c' = — \np, so that the condi- 
tion z = c/d = 1 reduces to hip ~ —2J/T. This agrees 
with the low-temperature asymptotics of the NL, whose 
exact form is e 2J / T = p/(l —p). It is thus plausible that 
the boundary of the Griffiths phase with divergent DoS 
in fact coincides with the NL as was conjectured in [l^ . 

The boundary at large disorder is set by the phase 
transition line into the paramagnetic phase, where the 
string mechanism ceases to be operative. Indeed, our 
numerics indicates pseudogap formation everywhere in 
the paramagnetic phase. Right at this part of the phase 
boundary (between the Nishimori point N and the per- 
colation critical point P c ), we find a divergent DoS with 
a non-universal exponent z, Fig. While we do not 
have at present an analytical theory for this behavior of 
DoS at criticality, an analogy with the random phase XY 
model, which has a phase diagram remarkably similar to 
Fig.Q with a low-T part of the phase diagram requiring a 
strong-disorder RG treatment, |24| might be instructive. 

In conclusion, we have reported the observation of a 
novel type of insulating phase for 2d electron systems 
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